Coexistence of Singlet and Triplet Attractive Channels in the Pairing Interactions 

Mediated by Antiferromagnetic Fluctuations 
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We propose a phase diagram of quasi-low-dimensional type II superconductors in parallel magnetic 
fields, when antiferromagnetic fluctuations contribute to the pairing interactions. We point out 
that pairing interactions mediated by antiferromagnetic fluctuations necessarily include both singlet 
channels and triplet channels as attractive interactions. Usually, a singlet pairing is favored at 
zero field, but a triplet pairing occurs at high fields where the singlet pairing is suppressed by the 
Pauli paramagnetic pair-breaking effect. As a result, the critical field increases divergently at low 
temperatures. A possible relation to experimental phase diagrams of a quasi-one-dimensional organic 
superconductor is briefly discussed. We also discuss a possibility that a triplet superconductivity is 
observed even at zero field. 

PACS numbers: 74.70.Kn, 74.20. Mn, 74.20.-z 



Pairing interactions mediated by a spin fluctua- 
tions have been discussed as a possible mechanism of 
anisotropic superconductivity (lj-|9|] in exotic supercon- 
ductors such as organic superconductors, heavy fermion 
superconductors, high-T c cuprate superconductors, and 
ruthcnatc superconductors. For example, in a ruthcn- 
ate superconductor Sr2Ru04, pairing interaction medi- 
ated by ferromagnetic fluctuations is a candidate for 
the mechanism of the superconductivity. The supercon- 
ductivity in this compound is considered to be due to 
triplet pairing from experimental results such as those 
obtained in Knight shift measurements |l0| ] and /iSR ex- 
periments p| . 

On the other hand, in the organics, cuprates, and 
heavy fermion superconductors, anisotropic singlet pair- 
ing intereactions mediated by antiferromagnetic fluctu- 
ations have been examined by many authors, because 
these compounds are in proximities to antiferromagnetic 
phases. However, in some of the quasi-one-dimensional 
organic superconductors and the heavy fermion super- 
conductors, triplet superconductivities have been sup- 
ported by some experimental results JlJ,[l3|]. In partic- 
ular, in some of the heavy fermion superconductors, the 
triplet pairing seems to be established. 

In a quasi-one-dimensional organic superconductor 
(TMTSF) 2 PF 6 , triplet pairing is supported by recent 
Knight shift measurements by Lee et al. |ll| , but it may 
appear to contradict the proximity to an antiferromag- 
netic phase. We will briefly discuss later that it does 
not necessarily contradict the mechanism of supercon- 
ductivity with pairing interactions mediated by antifer- 
romagnetic fluctuations. There is another contradiction 
between the thermal conductivity measurements [ [l2f and 
the NMR experiments @|5| in (TMTSF) 2 C10 4 . The 
former indicates a nodeless gap, while the latter suggests 
an existence of line nodes. However, this contradiction 
can be consistently explained also in terms of the pair- 



ing interactions mediated by antiferromagnetic fluctua- 
tions @. 

In the cuprate high-T c superconductors, experiments 
with 7r-junction suggest an anisotropic pairing which is 
conventionally called a d-wave pairing. Pairing interac- 
tions mediated by spin fluctuations have been discussed 
by many authors since the discovery of the high-T c super- 
conductors, for example, Miyake et al. M and the present 
author et al. 0. However, it is not established to what 
extent the spin fluctuations contribute to the pairing in- 
teractions. 

In this paper, we discuss a phase diagram of quasi- low- 
dimensional type II superconductors on a T-H plane. We 
consider a situation in which a singlet pairing is favored 
at zero field because of pairing interactions induced by 
the exchange of antiferromagnetic spin fluctuations. We 
restrict our discussion to parallel magnetic fields so that 
we can regard the orbital pair-breaking effect as weak. 

The pairing interactions mediated by spin fluctuations 
are written as 



(i) 



where i and j denote lattice sites and Xij is the expression 
of the spin susceptibility in real space. The dynamical 
effects and strong coupling effects can be taken into ac- 
count by developing a perturbation theory in more basic 
models such as Hubbard models However, since it 

is not essential for the purpose of this paper, we discuss 
pairing interactions of the form of eq. (1) for a while. 
Equation (1) can be rewritten as 
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(2) 
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where 



ipu (k) = c_ kT c kT 

V>i-i(k) = c_ ki c ki 

^io(k) = ^(c- kT c ki + c_ ki c kT ) 

V'oo(k) = -^(c_ kT c ki - c_ kJ c kT ). 



(3) 



When the antiferromagnetic fluctuations are strong, 
the factor \ij changes its sign depending on the distance 
between the sites i and j. In general, the factor Xij is pos- 
itive for the sites (i,j) on the same sublattice, while Xij 
is negative for the sites on the different sublattices, 
for the antiferromagnetic correlations. Here, it is easily 
verified from eq. (2) that the spin correlations between 
the different sublattices (xij < 0) favor singlet pairing, 
while they on the same sublattices (x%j > 0) favor triplet 
pairing. 

For example, in tight binding models near the half- 
filling, the factor \ij changes its sign alternately. For the 
nearest- neighbor sites Xij = Xi < 0, while for the 

next-nearest-neighbor sites Xij = X2 > 0. The former 
Xi < favor the singlet pairing, while the latter X2 > 
contribute to the next-nearest-neighbor triplet pairing. 

The triplet pairing interactions between the next- 
nearest-neighbor sites have not been considered to be 
important so far, because the nearest-neighbor singlet 
pairing interactions are usually much stronger than the 
next-nearest-neighbor triplet pairing interactions. How- 
ever, we could not ignore the triplet pairing interac- 
tions between electrons with parallel spins at high fields 
where the singlet pairing is suppressed by the Pauli pair- 
breaking effect. 

It might be considered that the contributions to the 
triplet pairing interactions from xi > may be com- 
pletely canceled by the larger negative contributions from 
X\ < 0. In order to illustrate that an attractive channel 
remains as a total, let us consider an example in which 
there are only two kinds of Xij with opposite signs, and 
write them xi arL d Xi- F° r example, in a one-dimensional 
tight binding model, if we truncate Xij &t the next- 
nearest-neighbor sites, we obtain xi < and xi > 
and |xi| > \Xz\- 

Then, the gap equations are written in a matrix form 



A 2 



= -V 



XiWn 
X2W21 



X1W12 
X2W22 



Ai 
A 2 



(4) 



where V = 3g/8 > for the singlet pairing and V = 
—g/8 < for the triplet pairing. The matrix elements 
W nm are defined by 



1 



tanh ■. 



^ m = -^ 7 „(k) = - 

k 



7m (k), 



(5) 



where 7n(k) are the momentum dependent factors to ex- 
pand the gap function. For example, when we consider 



one dimensional systems with interactions up to the next- 
nearest- neighbor sites, 7n(k) are defined by 



71 (fc^ = \Z2cosfcz 

72 (M = v / 2cos2fc a; 



for the singlet pairing, and 



Ji(k x ) = V^sinfca; 
l2(k x ) = v / 2sin2fc I 



(6) 



(7) 



for the triplet pairing. Extensions to two dimensional 
systems are straightforward. For example, j n are defined 
by 



71 (k) = cos k x — cos k y 

72 (k) = cos 2k x — cos 2k y 

for a d-wave-like pairing, while they are defined by 



71 (k) 

72 (k) 



y/2 sin k x 
2 sin k x cos fc„ 



(8) 



(9) 



for a p-wave-like pairing. The matrix elements W11 = 
W21 are not equal to zero for non-half-filling. 

Now, we show that the two eigen values of the matrix 
on the right hand side of eq. (4) have opposite signs in 
general. They are equal to 



1 



A± = - [(xiWii + X2W22) 



± 



\J {XiWn + X2W22) 2 - 4x1X2(^11! 



w 2 12 ) 



vW 2 2 

(10) 

Since X1X2 < and W11W22 — W 2 2 > 0, the second 
term in the square root is positive and thus the two eigen 
values A± have opposite signs. W11W22 
be proved as follows. First, we consider an average 



W12 > can 



1 x ^ tanh 
N 2^ o, )■ 



2e k 



(11) 



Then, since an inequality ((71 a; — 7 2 ) 2 ) > is satisfied 
for arbitrary real number x, we obtain 

(7i72) 2 -(7i 2 )<72 2 ) <0. 

Thus, WnW 2 2 - WI2 > is proved. 

The fact that the two eigen values have opposite signs 
means that both the singlet and triplet interactions have 
attractive channels. For antiferromagnetic fluctuations, 
the attractive triplet interaction is weaker than the at- 
tractive singlet interaction as we discussed above. Thus, 
usually, the singlet pairing should occur at zero field. 
However, at high fields where the singlet pairing is sup- 
pressed, the triplet pairing takes place. With an assump- 
tion that magnetic field of the order of fi e H ~ Ao does 
not change the feature of the spin fluctuations drastically, 
we could predict a phase diagram as shown in Fig.^ on a 
T-H plane. 
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FIG. 1. Schematic phase diagram on a temperature and 
magnetic field plane. The vertical solid line at high fields 
shows the transition temperature of a triplet superconductiv- 
ity of parallel spin pairing. The solid curve shows the critical 
field of singlet superconductivity, which is lowered into the 
dotted line below the critical field of the triplet supercon- 
ductivity at low temperatures. The dotted-broken line is the 
transition temperature of the triplet superconductivity when 
the singlet pairing is suppressed. 

Since the magnetic field is parallel to the highly con- 
ductive layers, the orbital pair-breaking effect is not ef- 
ficient. Thus, the transition temperature of the triplet 
superconductivity of parallel spin pairing is reduced only 
slightly by the magnetic field. Hence, the triplet super- 
conductivity appears at high fields. 

The upturn of the critical field (d 2 H c /dT 2 > 0) of sin- 
glet superconductivity at low temperatures in Figffl is due 
to an appearance of Fulde-Ferrell-Larkin-Ovchinnikov 
state (FFLO state or LOFF state) 0|Lf]. In low dimen- 
sional systems, the FFLO state is enhanced by Fermi- 
surface effects at low temperatures jO^C^. A reentrant 
transition to a triplet superconductivity might be ob- 
served at higher fields pi] ], but the transition temper- 
ature could recover only up to the zero field value at 
most. 

Our schematic phase diagram is consistent with 
the experimental results in (TMTSF) 2 C10 4 Q and 
(TMTSF) 2 PF 6 @. The observed critical field also 
shows an upturn at low temperatures (d 2 H c /dT 2 > 0), 
and tends to increase divergently in a low temperature 
limit (8H c /dT — > -co), although the observated transi- 
tions have widths and are rather ambiguous. 

In the FFLO state, the singlet and the triplet order 
parameters are mixed in general, since the finite center- 
of-mass momentum q of the pairs breaks a symmetry in 
the real space, and at the same time the magnetic field 
breaks the rotational symmetry in the spin space. Hence, 
the FFLO state includes a component of a triplet super- 
conductivity of antiparallel spin pairing. Matsuo et al. 



examined the mixing of the order parameters in an FFLO 
state, and found that the critical field is remarkably en- 
hanced |2^]. In particular, the tri-critical temperature of 
the normal, BCS, FFLO phase is enhanced. 

The lower critical fields of the FFLO state and the tri- 
critical point are not drawn in Fig.[j], since it is difficult 
to assert their locations. It is also difficult to predict the 
orders of the transitions at the lower critical fields. In 
(TMTSF)2C104, since first order transitions below the 
upper critical fields have not been observed, the FFLO 
state might be suppressed in this compound by impurity 
scatterings or some other reasons. In such a situation, the 
upturn of the critical field does not occur in our phase 
diagram. Such a suppression of the FFLO state might 
be a reason of the widths of the transitions at the upper 
critical fields in (TMTSF) 2 C10 4 . 

In our phase diagram, we should be careful about 
the choice of the ratio of the transition temperatures of 
the singlet pairing and the triplet pairing at zero field. 
The transition temperatures are extremely sensitive to 
the strengths of the pairing interactions, Coulomb re- 
pulsive interactions, impurity scatterings for anisotropic 
pairings, and so on. For example, if the triplet pairing 
interactions are too weak although they exist, triplet su- 
perconductivity at high fields is not observed in practice. 

On the other hand, the result of the recent Knight shift 
measurements in (TMTSF^PFg Jl3| might be explained 
by taking a larger ratio of the transition temperatures. In 
this compound, triplet superconductivity is supported by 
the Knight shift measurements in a strong magnetic field 
B = 2.38[T], and the superconducting transition tem- 
perature was still about 40% of the zero field transition 
temperature T c ^ = 1.18[K] in spite of the high magnetic 
field. If the superconductivity at zero field is due to a sin- 
glet pairing mediated by antiferromagnetic fluctuations, 
the strong magnetic field B — 2.38 [T] may favor a triplet 
pairing, since the pairing interactions include attractive 
triplet channels. 

We have illustrated the existence of an attractive 
triplet channel in the interactions mediated by antifer- 
romagnetic fluctuations in restricted cases. However, it 
is plausible that this feature is general. For example, the 
interactions between distant sites were calculated in a 
perturbation theory based on the quasi-one-dimensional 
Hubbard model, and it was found that triplet pairing 
interactions can be attractive H . 

In a 1/ 4-filled Hubbard model with electron dispersion 



ek = — 2t cos k x a — 2t' cos k y b — /i, 



(12) 



a perturbation theory within an RPA gives pairing in- 
teractions V(Rij — ma + nb) = V mn shown in Table | 
and [[J, where a and b are lattice vectors. In these tables, 
we find oscillations of the interactions in the real space. 
For triplet pairing, for example, V\q < and V40 < 
are attractive interactions. The former favors the intra- 
chain nearest-neighbor component of the gap function 
A10 sm{k x a), which does not have a node on the open 
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Fermi-surfaces. On the other hand, the latter favors the 
gap function of the form A40 sin(4fc x a), which has line 
nodes on the Fermi-surface near k = (±7r/4a, ±ir/2b) for 
1/4-filling. In practice, a superposition up to more dis- 
tant components gives the gap function localized around 
the Fermi-surfaces with a width related to the spin fluc- 
tuations. 

TABLE I. Singlet pairing interactions Vnm/t in a 
quasi-one-dimensional Hubbard model with t' — 0.150t and 



U w 0.1481 





n = 


n= 1 


n = 2 


m = 


2.769 


-0.067 


0.020 


1 


0.181 


0.006 


-0.002 


2 


-0.465 


0.087 


-0.031 


3 


-0.084 


-0.004 


0.004 


4 


0.252 


-0.102 


0.036 


5 


0.056 


-0.002 


-0.005 



TABLE II. Triplet pairing interactions V nm /t in a 
quasi-one-dimensional Hubbard model with t' = 0.150* and 



U w 0.148*. 





n = 


n= 1 


n = 2 


m = 


0.768 


0.024 


-0.013 


1 


-0.078 


-0.002 


0.001 


2 


0.187 


-0.031 


0.010 


3 


0.032 


0.002 


-0.001 


4 


-0.100 


0.036 


-0.012 


5 


-0.020 


0.001 


0.002 



In conclusion, we point out that the pairing interac- 
tions mediated by antiferromagnetic fluctuations neces- 
sarily accompany attractive channels of triplet pairing 
interactions. Therefore, in quasi-low-dimcnsional type II 
superconductors in magnetic fields parallel to the highly 
conductive layers, if the origin of the pairing interac- 
tions is the exchange of antiferromagnetic fluctuations, 
a triplet pairing occurs at high fields where the singlet 
pairing is suppressed by the Pauli pair- breaking effect. 
We propose a schematic phase diagram which may ex- 
plain the experimental observation in (TMTSF) 2 C104. 

We have discussed quasi-low-dimensional systems in a 
strong parallel magnetic field. However, if some extra 
mechanism suppresses the singlet pairing more than the 
triplet pairing, antiferromagnetic fluctuations can give 
rise to a triplet pairing superconductivity even at zero 
field, or even in three dimensional systems. For example, 
intersite Coulomb repulsions may suppress the pairing 
between electrons on nearest- and next-nearest-neighbor 
sites. Also, the suppression of anisotropic supercon- 
ductivity due to impurity scatterings would be differ- 
ent depending on the pairing anisotropy. The quasi-one- 
dimensional organic superconductors (TMTSF^X have 
open Fermi-surfaces. Thus, the triplet pairing of the form 
A(k) ~ sm(fc x a), which was discussed by Belin et al. [ p"2| , 
does not have a node on the Fermi-surface. Such a node- 
less gap function might have some advantage in com- 
parison to the other anisotropic gap functions with line 
nodes. These possibilities will be examined in a separate 
paper. 

This work was supported by a grant for CREST from 
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